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B Meson Decay Constants Using NRQCD 

Junko Shigemitsu a * 

a Physics Department, The Ohio State University, Columbus, Ohio 43210, USA. 

Recent results for B meson decay constants with NRQCD b-quarks and clover light quarks are discussed. 
Perturbative matching factors through O(afM) are now available and incorporated into the analyses. An O(aa) 
improvement term to the heavy-light axial current is identified and included. The slope of fpsVMps versus 
1/Mps is significantly reduced by these corrections. 



1. Introduction 

The heavy-light pseudoscalar meson decay con- 
stant, fps, is defined through the matrix element 
of the heavy-light axial vector current between 
the pseudoscalar meson state and the hadronic 
vacuum. In Euclidean space one has, 



[0\A M \PS)=p li f PS . 



(1) 



fu 3 has been measured experimentally, giving 
values consistent with lattice predictions. /g, on 
the other hand is unlikely to be measured directly 
any time soon via leptonic B decays, B — ► I 77. 
This fact coupled with the importance of fp in 
analyses of BqBq mixing phenomena, makes an 
accurate lattice determination of fp particularly 
relevant H. 

This talk describes work being carried out by 
A. Ali Khan, T. Bhattacharya, S. Collins, C. 
Davies, R. Gupta, C. Morningstar, U. Heller, J. 
Sloan and myself, on heavy meson decay con- 
stants with NRQCD b-quarks. Earlier works ex- 
ist by C. Davies {§, Draper & McNeile f| and 
S. Hashimoto Q]. Dr. Onogi discusses the Hi- 
roshima group's results on NRQCD decay con- 
stants in a separate talk ||. 

New developments since Lattice '96 include, 

• Completion of the one-loop matching cal- 
culation between lattice NRQCD and con- 
tinuum full QCD axial currents through 
0(a/M). 
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• Inclusion of an O(aet) discretization correc- 
tion to the local heavy- light axial current. 

• Considerable decrease in the slope of 
fpsVMps versus 1/Mps once the dis- 
cretization correction and the matching Z- 
factors have been included. This makes 
the connection between the static limit and 
the physical b-quark region much smoother 
than previously thought. 

• Higher statistics on quenched calculations 
at (3 = 6.0 that include 1/M 2 corrections to 
the action and currents at tree level. This 
latter data is still being analyzed and I will 
only show older data in this talk with cor- 
rections through 0(1/M). 

2. The Action and Current Operators 

We use the NRQCD action to simulate b-quarks 
It is given by (in continuum notation), 



i> (D t + H + SH) i\>, 



(2) 



We have dropped the rest mass term and Hq is 
the nonrelativistic kinetic energy operator, 

£>(2) 



H 



2M ' 



and through 0(Aqcd {^qcd/M)) one has 
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In our more recent simulations we have added the 
terms 



8(M y 



rcr-(DxE-ExD) 
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which are of 0(Aq C d (Aq C d/M) 2 ) and 

P (2) ) 2 f5) 

8(M ) 3 1 j 

which we believe is the domi- 
nant 0(Aqcd (Aqcd/M) 3 ) contribution. These 
latter calculations also include discretization cor- 
rections to the lattice laplacian and the lattice 
time-derivative 0. All operators in the NRQCD 
action and in the heavy-light currents discussed 
below, are tadpole improved (£/—>£/ /uq with 
uq = fourth root of the average plaquette). Tree 
level coefficients are used in the action. 

For light quarks we have data using the Wilson 
and the Csw = 1 clover quark actions. The bulk 
of our results, however, come from the tadpole 
improved clover action. Gauge fields were created 
using the Wilson plaquette action. The quenched 
configurations are at /? = 6.0; one set provided by 
the UKQCD collaboration and another newer set 
created by us on the LANL and NCSA CM5 ma- 
chines. The dynamical calculations were carried 
out on the HEMCGC nj — 2 staggered configu- 
rations at (3 = 5.6. 

Heavy-light currents in full QCD have the form 
qTh. The four component Dirac spinor for the 
heavy quark, h, is related to the two component 
NRQCD heavy quark (heavy anti-quark) fields, 
■0 (ip), via an inverse Foldy-Wouthuysen trans- 
formation. 



h 



u 



FW 



(G) 



TJ- 1 

U FW 



1 H ■ D) 

2M u ; 

5(2 



8M- 



{D (Z > + gH ■ B — 2iat ■ E) 

+ 0(l/M 3 ) (7) 



a = 707 and S = diag{cr, or). All our expressions 
have been converted to Euclidean space with her- 
mitian 7-matrices. 

Through 0(1/M) one has at tree-level, 
j = j(o) +J (D 



qTQ - 2i??r(7 • D)Q 



(8) 



In the next section we discuss what happens at 
one-loop. 

3. Perturbative Matching 

In order to match between lattice currents used 
in our simulations and those of continuum QCD, 
we consider the process in which a heavy quark of 
momentum p is scattered by the heavylight cur- 
rent into a light quark of momentum p' . For a 
one-loop matching we need to go through the fol- 
lowing steps. 

1. Carry out the one- loop calculation for the 
above process in full QCD. 

2. Expand the amplitude in terms of p/M , 
p'/M etc. 

3. Identify operators in the effective theory 
(NRQCD) that would reproduce these 1/M 
corrections. 

4. Carry out a one-loop mixing matrix calcu- 
lation in the effective theory. 

For HQET analogous calculations have been done 
by Eichten & Hill §, Golden & Hill §,and Neu- 
bert [[To) . Colin Morningstar and I have now com- 
pleted the matching calculation, between lattice 
NRQCD and full continuum QCD, for the time 
component of the axial vector current through 
0(a/M) |~L|. 



Continuum calculation 

A one-loop calculation in full QCD finds that 
(q(p') I Aq I h(p) )qcd involves terms propor- 
tional to, 



u q (p')l5loUh{p) 



M 



U q (p'h5Uh(p) 



^f- U g {p')jBlOU h (p) , ll U q {p')~f 5 U h {p) 

P ■ P' Po 
M 2 M 



W)l^h{p) + 0(1/M 2 ) (9) 



We set the light quark mass equal to zero, i.e. we 
ignore terms such as m q /M . 
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The above terms can be reproduced via three op- 

(i) - (i) 

erators in the effective theory, J\' — qO A Q , 
i = 0,1,2 , after making use of equations of motion 
for the light quark. 



J 



(o) _ 



(i) 



(2) 



= 975 7o Q 

9757o(7 D)Q 
(Dq ■ 7)7570*3 



(10) 



ap- 



[q(p')\A \h(p) 



'QCD 



can be written as, 



QCD 



A 1 \ 

A I ren 



i=0 
2 



YsVffaip'^Ugip) 



(11) 



a local heavy-light current, and J AL and J A \ 
have D — > {lattice symmetric covariant deriva- 
tive}. These lattice operators are only defined 
up to improvement terms. In fact, in the next 
section we will argue that a consistent one-loop 
calculation with clover (i.e. 0(a) improved) light 
quarks, requires that an O(aa) correction term 



be added to the local J 



(o) 

A,L ' 



The lattice matrix elements ( J 



(0 » 

A,L 1 



can be re- 



lated to the ( J 
a mixing matrix Z. 



on the RHS of eq.(hlf) via 



J(°) and JW coincide with eq.(|), J*- 2 

pears only at one-loop. So, through one-loop ( J/ll ) = ^ Z,j ( J 



C?) 



(13) 



From (0) and (|l|) one extracts the matching 
relation between matrix elements in full QCD and 
those evaluated in lattice simulations, 



CO 



! ^0 )q C D - ^2 vfZij 1 ( J A, 

id 



T<J) > 

'a,l , 



(14) 



6$ are momentum space representations of the 
operators between quark fields in eq. (|ic|) . We 
have calculated the one-loop coefficients, r]f, us- 
ing dimensional regularization (with totally anti- 
commuting 75) in the MS scheme. A gluon mass, 
A, is introduced as an infrared regulator. We find, 
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(12) 



The infrared divergent terms will cancel when we 
match to a lattice regularized one-loop calcula- 
tion. 

Lattice calculation 

The lattice currents used in our simulations are 
discretized versions of the operators in eq. (|l^) . In 
the absence of any improvement, J AL becomes 



C = ff Tr 1 is the vector of matching coefficients 
for this problem. 

We have calculated the one-loop contributions to 
Zij in lattice perturbation theory for lattice ac- 
tions including H and SH of eq. ([!]). Equation 
( PH ) can be expanded out as, 

(A) )qcd = 

{1 + a [B - \(C q + C Q ) - Coo - Go]) {Jfl) 
+ (l + a[B 1 - + C Q ) - C01 - C11] ) (J { A) 

+ a[B 2 - C02 - C12 ] (J ( aI) + 0(a 2 ) (15) 

The BiS come from the nf's in ( |l2"| ) and C q 
and Cq are the light and heavy quark lattice 
wave function renormalizations respectively. The 
are the one-loop vertex correction contribu- 
tions to Zij. Before proceeding with combining 
perturbative numbers with simulation results for 
the matrix elements ( J 1 1 \ ) , we must discuss im- 



provement of the local current J AL 
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4. An O(ao) Correction to the Heavy- 
Light Axial Current 

There has been a lot of work recently on im- 
proving quark bilinear operators such as the vec- 
tor and axial currents. For the light quark sector, 



the DESY 12 group has shown that an improved 
axial vector current takes on the form (suppress- 
ing isospin indices), 



A 



A 



toe 



c A aDnP 



(16) 



with P the pseudoscalar density and the sym- 
metric lattice derivative. In perturbation theory, 
ca starts out at 0(a) flj,[l3| ■ Less is known about 
improvement of the local heavy-light or static- 
light axial currents. Borrelli & Pittori jhij have 
shown that through one-loop, static-light bilin- 
ears using clover light quarks are free of 0(a) 
and 0(a alog{a)) terms. They did not consider 
0(a a) terms. 

Our one- loop calculation of the mixing matrix Z^j 
shows that there is an 0(a a) lattice artifact term 
whose effects can be removed by improving J A °\. 
The improved heavylight lattice axial current ac- 
quires a correction term that is the precise ana- 
logue of the correction term in eq. (|l6|) . The ar- 
gument goes as follows. 



In order to calculate Zn 



a C02 one needs to 



start from ( J A \ ) ano - project out terms propor- 

(21 l 

tional to ( J\ ) ren = ^ u q (-ip' ■ 7)7570 Uq. 
After calculating (02 m this way, we find that £02 
has a term that grows with a M, the dimension- 
less heavy quark mass. So, 

a C02 ( J A 2) )ren = « (dM A X + A 2 ) ( ) ren 

has an 0(ao) term with the structure, 

aau q (-ip' ■ 7)7570 U Q . 
One can interpret this term as coming from a lat- 
tice operator, 



j ax = a (Dg • 7)7570*2 
and what we are finding is, 



(17) 



a C02 ( J A ^ )re-< 
_ /-true 

= "C 02 



J 



(2)- 



+ aCo 3 (J (3) 



A,L l Ten 



(18) 



and from eq.(0 



£ ■ 

j=0,l,2 



J 



U)- 



+ 01 C03 ( J A L )ran 



(19) 



with Zq2 = a Co2 Ue ■ K one wants to have the same 



set of matrix elements on the RHS of eq. (|19|) as in 
the continuum theory, the last term must be re- 
moved. This is easily accomplished by improving 



J 



(o) 

A,L- 



t(0) 
J A,L 



J A,L 



- J A,L - a( *03-J A ,L 



and one now has 



J 



.1.1 z2 ■ 

.7=0,1,2 



} 



J 



(./) 



(20) 



(21) 



Going through the same steps that led from ( |l3| ) 
to ([u]) and then to (|l5|), one ends up with an 
equation very similar to eq.(|i"5|) with 



'^°L>->(4 0) />and Co2-a7 e - 



A,L 



(2) (3) 

Although J A ' L and J A ' L are proportional to each 

they play very different 



other, J { A \ = 2,iM J 

r(2) 



(2) 

A,L ' 



roles. J A L is the lattice version of a current op- 
erator that exists in the continuum theory. It is a 
l/M correction to the static heavy-light axial cur- 

(3) 

rent and is absent in the static theory. J A L , on 
the other hand, has no continuum counter part. 
It survives into the lattice static theory. So the 
static limit of eq.([l5|) becomes, 

(Aq )qcd 

= (l + a[B - \(C q + C Q ) - Coo]) (Jj^/W 
= (1 + a [Bo - \(C q + C Q ) - Coo]) (J A %)stat 



" Co3(J ( A)stat + 0{a 2 ) 



(22) 



We will see in the next section that the last cor- 
rection term in (^2|) significantly reduces the value 
of fpsVMps m the static theory. 

(3) 

It is easy to see that the improvement term J A T 
is of the same form as the second term in eq.(|16|). 
If one defines a heavylight pseudoscalar density 
Phl = 975*2, then 
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(3) 

gives the operator J\ ' L of eq. ( |17[ ) 



aaD P H L = aa(D q)^5Q + 0(aa/M) 

w -aa (D 9 7o)757oQ (23) 

In the above we ignore the term coming from the 
time-derivative acting on the heavy quark field 
Q, since equations of motion make that into an 
0(a a/M) term and we neglect such contributions 
together with contributions of 0(a/M 2 ). Apply- 
ing light quark equations of motion to the last ex- 
pression in ( |25| ) 
multiplied by a 



5. Some Quenched fs Results 

An analysis of NRQCD heavy meson decay 
constants including one-loop matching factors, 
has been carried out for the first time in Ref . jl5| . 
The data was obtained on 16 3 x 48 quenched con- 
figurations at (3 = 6.0. Both the gauge configura- 
tions and light propagators were generously pro- 
vided by the UKQCD collaboration. Tadpole im- 
proved clover light propagators were used. The 
NRQCD action included the SH of eq.(f§) but 
no 1/M 2 terms. Two light n values around the 
strange quark mass and four NRQCD bare heavy 
quark masses were used. We also have results for 
static heavy quarks. 

For the perturbation theory, we use ay of Lep- 
agefcMackenzie |L6j . The q* for this matching 
calculation has not been calculated yet. So, we 
present results for both q* = 1/a and q* = ir/a. 
The log{aM) terms that appear in the matching 
coefficients after cancellation of logarithmic IR di- 
vergences, are set to a constant value log(amb) for 
all bare heavy quark masses, nib = 4.1(l)GeV is 
the b-quark MS mass at a scale equal to its value 
|l7| . Alternatively we could have used the one- 
loop renormalization group improved expression 

with an overall factor of (a(mb) / ct^a -1 )) 2 ^° for 
the two matching coefficients Co and C\ . The dif- 
ference between the two approaches is very slight 
and at the B-meson numerically undetectable. 
Since we presently have results at only one lattice 
spacing, we have used the simpler log(timb) pre- 
scription. Details of the simulations and fitting 
procedures are given in Ref. p5[ . Figure 1. shows 
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Figure 1. Heavy-Light Decay Constants at n s . 
Circles : tree-level ; Diamonds : one-loop match- 
ing with aq* = ir ; Squares : one-loop matching 
with aq* = 1. We demonstrate the effect of the 
correction term to L , by showing the one- loop 
matched result in the static limit without the last 
term in cq.(p2|) for aq* = n (upper burst) and 
aq* = 1 (lower burst). The physical B s meson is 
just below 1/aMps = 0.4. 



a 3 ^ 2 fpsVMps versus l/aM PS for k — K s . We 
show the tree-level results and the one-loop re- 
sults for the two q*'s. For the static limit we also 



indicate the one-loop values without the ( J 



(3) 

A,L , 



correction term. The physical B s meson is just 
below 1/aMps = 0.4 on the figure. 

There are several interesting things to note in Fig- 
ure 1. 

1. After including the one- loop matching fac- 
tors and improvement of J^\, the slope of 
a 3 ^ 2 fps^Mps versus 1/aMps as one leaves the 



G 



static limit decreases considerably. For q* = 1/a 
it is consistent with zero. For q* — ir/a the rela- 
tive slope ( the slope divided by (fpsVMps) tat 
) is ~ — lGeV. Unfortunately a precise value of 
q* is required in order to make a more quantita- 
tive statement. Ref. []l8| quotes q* = 2.18/a for 
the static theory. 

2. Around the physical B meson region the one- 
loop corrections are a 9 - 13% effect depending on 
q* at k s . The difference due to q* of 1/a or 7r/a 
is a ^5% effect on the final answer. This is one 
measure of the uncertainty in fg s coming from 
higher orders in the matching calculation. 
The static limit is much more sensitive to q* and 
the one-loop corrections are large (25 - 40% at 
k s and 10 - 30% at re c ), with a third of the shift 
coming from the O(aa) correction to J^ L - 

We use a scale of a -1 = 2.0(2)GeV to convert to 
physical units. In a quenched calculation, differ- 
ent observables can lead to different a _1 's. We 
have used an a -1 consistent with quenched light 
quark calculations (M p , etc.) and allow for a 
large error of 10%. This feeds back into one of the 
dominant systematic errors in our final estimate 
for fs- We find at k s , 



Sb. = 



0.198 (8)(30)(17) GeV 
0.209 (8) (32) (17) GeV 



and after extrapolating to k c 



Sb = 



0.174 (28)(26)(16) GeV 
0.183 (32)(28)(16) GeV 



1/a 
ir/a 



1/a 
7r / 'a 



Errors in the first brackets correspond to statis- 
tical plus fitting plus k extrapolation errors. In 
the second brackets we give a -1 systematic errors 
and the third brackets summarize our estimates 
for higher order perturbative and 1/M 2 correc- 
tions. We now have more recent simulations with 
1/M 2 terms included at tree level both in the 
NRQCD action and the currents (see next sec- 
tion on Work in Progress). Their contributions 
are at the 3 - 4% level. The Hiroshima group 
reports similar findings for the 1/M 2 current cor- 
rections ||. We have not included any estimate 
for continuum extrapolation corrections, since at 
the moment we only have results at a single value 



of 8. Studies with Wilson and/or static fermions 
|l9| ||l| have shown noticeable lattice spacing de- 
pendence around 3 = 6.0. It will be interesting 
to do a thorough scaling study with clover light 
fermions and an improved J^l- 



6. Work in Progress and Future Plans 

Two sets of NRQCD heavy meson decay con- 
stant data are currently being analyzed by my 
collaborators. 

1. The Glasgow-LANL-OSU-Kcntucky (GLOK) 
collaboration has results on 16 3 x 48 quenched 
configurations at 8 = 6.0 §§]. Both thc NRQCD 
action and the current operators include 1/M 2 
corrections at tree level. We also use a better 
time evolution equation to obtain heavy propa- 
gators ( which follows from a slightly modified 
lattice NRQCD action) than in the simulations 
of the previous section. The new evolution equa- 
tion eliminates a residual 0(aA(A/M)) error in 
amplitudes. The light quark action employed is 
the tadpole-improved clover action. We use five k 
values and six heavy quark masses. Perturbative 
matching calculations for the new lattice action 
have now also been completed JTT| . Some prelimi- 
nary tree-level results were presented at St. Louis 



2. The SCRI-Glasgow-OSU (SCO) collabora- 
tion has results on the HEMCGC n/ = 2, 
arridyn = 0.01 dynamical staggered configura- 
tions at 8 = 5.6. Both Wilson and tadpole- 
improved clover light quarks have been combined 
with NRQCD heavy quarks. The NRQCD ac- 
tion and the heavy-light currents were corrected 
through 0(1/M). Three k values and eleven 
heavy quark mass values were used. The Wil- 
son light fermion results are published in Ref. [M . 
Tree level analyses of the clover light quark data 
have been reported on in |2f|. We are now com- 
pleting the analysis including one-loop matching 

Several matching coefficient projects are also on 
the agenda. The one-loop matching calculation 
for the heavy-light vector current is well under- 
way and g*'s for both axial and vector currents 
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are high on our list. We have also started studies 
of nonperturbative renormalization of NRQCD 
operators. 

Finally, quenched simulations at = 5.7 have be- 
gun and in the future we plan to go onto (3 = 6.2. 

7. Summary 

The NRQCD approach to B meson decays is 
looking very promising. The l/M expansion is 
working well and appears to be under good con- 
trol. The first one-loop perturbative matching 
calculations have been completed and incorpo- 
rated into our analyses. Uncertainties due to 
higher orders in perturbative matching are es- 
timated to be at the ~5% level around the B 
meson. Many more results should be forthcom- 
ing soon, including the use of an improved heavy 
quark time evolution equation, and studies of 
l/M 2 corrections, unquenching and scaling. 
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